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Abstract: We consider neutrino oscillations in non- uniform matter in a quantum field 
theoretic (QFT) approach, in which neutrino production, propagation and detection are 
considered as a single process. We find the conditions under which the oscillation prob- 
ability can be sensibly defined and demonstrate how the properly normalized oscillation 
probability can be obtained in the QFT framework. We derive the evolution equation for 
the oscillation amplitude and discuss the conditions under which it reduces to the standard 
Schrodinger-like evolution equation. It is shown that, contrary to the common usage, the 
Schrodinger-like evolution equation is not applicable in certain cases, such as oscillations 
of neutrinos produced in decays of free pions provided that sterile neutrinos with Am 2 > 1 
eV 2 exist. 
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1 Introduction 

In experiments with solar, atmospheric and supernova neutrinos, the neutrinos propagate 
significant distances in matter before reaching detectors; the same will also be true for the 
proposed very long baseline accelerator neutrino experiments. Matter can affect neutrino 
oscillations drastically, leading, in particular, to their resonance enhancement through the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect [1, 2] or through the parametric enhancement 
mechanism [3]. It is thus extremely important to put the analyses of neutrino oscillations 
in matter on a solid theoretical basis. 

The standard approach to neutrino oscillations in matter, pioneered by Wolfenstein 
[1], is as follows. Mass eigenstate neutrinos fj, composing a given flavour neutrino state 
u a , are assumed to have the same momentum (and, due to their different mass, different 
energies). In the case of oscillations of relativistic neutrinos in vacuum, the evolution of the 
transition amplitude in the flavour basis is then described by the Schrddinger equation 1 
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is the neutrino mass matrix in the mass eigenstate basis, and U is the leptonic mixing ma- 
trix which relates the navour-eigenstate neutrino state vectors \v a ) {a = e, /i, r, . . . ) with 
the mass-eigenstate ones \vj) (i = 1,2,3, ...): 



Since the first term in the square brackets in the expression for Hq in (1.1) is proportional to 
the unit matrix, it leads to a common phase shift of all neutrino flavour states and therefore 
does not affect the oscillation probabilities. Thus, it can be omitted from Hq. The matter 
effects on neutrino oscillations can be incorporated by replacing the free Hamiltonian Hq 
by the effective Hamiltonian of neutrino propagation in matter according to Hq — )■ H = 
Hq+V. Here V is the matrix of matter- induced neutrino potentials due to coherent forward 
scattering of neutrinos on matter constituents; in the absence of background neutrinos it is 
diagonal in the flavour-eigenstate basis. Taking into account that for pointlike relativistic 
neutrinos the distance x they propagate over the time t satisfies x ~ t and that to leading 
order in small neutrino masses p ~ E where E is the average neutrino energy, one arrives 
at the following equation describing neutrino flavour evolution in matter: 



Equation (1.3) is employed in virtually all studies of neutrino oscillations in non- 
uniform matter. However, its derivation presented above was based on heuristic consid- 
erations and it certainly needs to be put on a more solid ground. Attempts at deriving 
eq. (1.3) within the relativistic quantum mechanics and quantum field theory (QFT) frame- 
works have been made in a number of papers. In [4] evolution of Dirac neutrinos in matter 
was described by making use of a Dirac equation with matter-induced potential. It was 
demonstrated that the neutrino wave function satisfies eq. (1.3) provided that matter den- 
sity varies little over the distances of order of the neutrino de Broglie wavelength. The 
Dirac equation was also employed in [5] , though only the case of matter of constant density 
was considered there. In ref. [6] the Dirac equation was used to describe the evolution of 
Dirac and Majorana neutrinos in matter, but again only in the case of constant-density 
matter. Besides the already mentioned ref. [4], the Dirac equation has been employed for 
describing neutrino oscillations in non-uniform matter in refs. [7, 8], both in the Dirac [7, 8] 
and Majorana [8] neutrino cases. In none of these papers, however, neutrino production 
and detection processes were included in the description of neutrino oscillations. The most 
advanced study of neutrino evolution in matter of varying density was carried out in [9] 
in the QFT framework. In that paper the treatment included the neutrino production 
and detection processes, and it was also demonstrated how the correctly normalized os- 
cillation probability can be obtained. The employed normalization procedure was rather 
cumbersome, though. 

The main goal of refs. [4-9] was to derive the evolution equation (1.3) from relativistic 
quantum mechanics or QFT. However, the conditions under which this equation is valid 
were not fully studied in those papers. Furthermore, the question of how neutrino os- 
cillations can be described in the situations when these conditions are not satisfied (and 
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consequently eq. (1.3) cannot be used) was not addressed. In addition, no discussion of 
neutrino production and detection coherence and of their effect on neutrino oscillations in 
matter was given. 

In the present paper we consider neutrino oscillations in non-uniform matter in the 
framework of QFT. In this approach neutrino production, propagation, and detection are 
treated as a single process, described by a Feynman diagram with the neutrino in the in- 
termediate state (such as the one in fig. 1). We discuss the conditions under which the 
oscillation probability can be extracted from the rate of the overall neutrino production- 
propagation-detection process and demonstrate that this probability is automatically cor- 
rectly normalized and satisfies the unitarity constraints. We also identify the conditions 
under which the amplitude of neutrino flavour transition can be found as a solution of 
eq. (1.3). One of our main results is that eq. (1.3) is not applicable when neutrino pro- 
duction and/or detection coherence is violated. We discuss the situations when this can 
happen and consider an alternative way of describing neutrino oscillations in those cases. 

Our treatment of neutrino oscillations in non-uniform matter closely parallels the treat- 
ment of neutrino oscillations in vacuum performed in ref. [10], but differs from the latter 
in a number of important aspects. The differences are mostly related to the properties of 
the neutrino propagator in non-uniform matter, which deviate significantly from those of 
the vacuum neutrino propagator and do not allow using some techniques that were applied 
to the neutrino propagator in vacuum. To make it easier to follow our treatment, let us 
briefly outline our main steps. 

• We consider the neutrino production, propagation and detection process described 
by the Feynman diagram of fig. 1. The external legs in this diagram correspond to 
the particles that accompany neutrino production and detection. These are either 
propagating particles or bound states, which are described by the suitable state 
vectors. The intermediate neutrino state is described by a propagator, which is found 
as a solution of the corresponding Dirac equation with matter-induced potential for 
neutrinos V(x). 

• Since this potential depends on the coordinate x, the system is not translationally 
invariant and the neutrino momentum is not conserved. As a result, the neutrino 
propagator in the momentum space depends on two momenta, p and p', rather than 
one. The amplitude of the overall process can be written as the integral over these 
two momenta, with the production and detection amplitudes <&p(p) and $£>(p') 
multiplying the momentum-space neutrino propagator S(E;p',p) in the integrand. 

• For the intervals of momenta over which the propagator varies significantly, the neu- 
trino production and detection amplitudes $p(p) and ^^(p') change very little. This 
allows one to greatly simplify the expression for the amplitude of the process. 

• We calculate the rate of the overall neutrino production-propagation-detection pro- 
cess and identify the conditions under which this rate factorizes into the neutrino pro- 
duction rate dT P a° d / ' dE , propagation (oscillation) probability P a p and the detection 
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Figure 1. Feynman diagram describing neutrino production, propagation and detection as a single 
process. 

cross section a p. When these conditions are satisfied, the oscillation probability can 
be extracted from the rate of the overall process by dividing the latter by dT^ od / dE 



• We reconstruct the oscillation amplitude from the expression for P a p and derive the 
equation that it obeys. We identify the conditions under which this equation coincides 
with eq. (1.3) and also discuss the situations when these conditions are violated and 
eq. (1.3) is not applicable. 

The described above program is realized in sections 2-4 of the paper; in section 5 
we summarize and discuss the obtained results. To make the paper self-contained, in 
Appendix A we briefly review the derivation of the neutrino propagator in non-uniform 
matter performed in [9], whereas in Appendix B we give a compendium of expressions for 
neutrino potentials in matter. Appendices C and D contain derivation of some results used 
in sections 2.2 and 4. 

2 The transition amplitude 
2.1 General formalism 

Consider the process of neutrino production, propagation and detection described by the 
Feynman diagram of fig. 1. We shall be assuming that the neutrino production process 
involves one initial state and one final state particle (besides the neutrino). Likewise, the 
detection process will also be assumed to involve only one particle besides the neutrino in 
the initial state and one particle in the final state. The generalization to the case of an 
arbitrary number of particles involved in the neutrino production and detection processes 
is straightforward and would just make the formalism more cumbersome without providing 
further physical insight. 2 

Let us first discuss the state vectors of the particles accompanying neutrino production 
and detection (external particles). In quantum theory, one-particle states of particles of 

2 As only one particle is assumed to be in the initial state of the production process, it must be unstable. 
This will be of no importance for us here, though. 
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type A can be written as 

\A) = J[dp]f A (p,P)\A,p), (2.1) 

where \A, p) is the one-particle momentum eigenstate corresponding to momentum p and 
energy Ea(p), and /^(p, P) is the momentum distribution function with the mean mo- 
mentum P. In eq. (2.1) we use the shorthand notation 

[dp] = ^ . (2.2) 

(2tt)3 v /2^(p7 

For particles with spin, the states \A) and |A, p) depend also on a spin variable, which we 
suppress to simplify the notation. 

Throughout this paper we will be using the normalization conventions of ref. [13] and 
the notation Pl,r = (1 =F 7s)/2. We choose the Lorentz invariant normalization condition 
for the plane wave states \A, p): 

(Ap'lAp) = 2£ A (p)(27r) 3 <5( 3 )(p-p / ). (2.3) 

The standard normalization of the states (-A|-A) = 1 then implies 

' I/a( P )| 2 = I- (2.4) 



(2tt) 3 

The states describing the external particles in fig. 1 can be represented in the form 
(2.1). For the initial and final states at neutrino production we write 

\ p i) = J [M fPiil, Q) \ p i, q) , \P f ) = J[dk]f Pf (k,K)\P f ,k), (2.5) 

and similarly for the states participating in neutrino detection: 

|A) = yi^/«(q',Q')lA,q'>, \D f ) = J[dk'}f Df (k',K')\D f ,k'). (2.6) 

We assume these states to obey the normalization condition (2.4). Some (or all) of the 
mean momenta of the external particles Q, K, Q' and K' may vanish, i.e. the states in 
eqs. (2.5) and (2.6) can describe bound states at rest as well as wave packets. 

The amplitude of the neutrino production - propagation - detection process is given 
by the matrix element 

iAt. = <r l D lf t~ P [-if**HM]-l\*Dt, (2.7) 

where T is the time ordering operator and Hi(x) is the charged-current weak interaction 
Hamiltonian. From this equation it is easy to calculate the transition amplitude in the 
lowest nontrivial order in T~Li using the standard QFT methods. The resulting expression 
corresponds to the Feynman diagram of fig. 1 and can be written as 

Uf, a = J[dq] fp^Q) J[dk}f* Pf (k,K) 

W] /«(q', QO / [dk'\ fhfQt', K') iA p ^-(q, k; q', k') . (2.8) 
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Here the quantity -A^'(q, k; q', k') is the amplitude of the process with plane- wave external 
states: 

iA^{q,k-q'X) = f d A Xl f d*x 2 M D tf,kf)e- i «- k 'X x *-*ri 

x S Pa (x 2 , Xl )M P (q, k) e -*(«-*)(*i-^) . (2.9) 

In this equation x\ and x 2 are the 4-coordinates of the neutrino production and detection 
points. The choice of the 4-coordinate dependent phase factors corresponds to the assump- 
tion that the peaks of the wave packets of particles involved in the production process 
are all located at xi = xp at the time t\ = tp, whereas for the detection process the 
corresponding peaks are all situated at x 2 = x^> at the time t 2 = to (this assumption can 
be relaxed, see section 6.2 of ref. [10]). The quantities Mp(q,k) and Mp>(q',k') are the 
plane- wave amplitudes of the processes Pi — > Pf + v a and Di + vp — > Df, respectively, with 
the neutrino spinors excluded. They are related to the full plane- wave neutrino production 
and detection amplitudes Mp(q,k) and Mp)(q',k') through 

M P (q, k) = ^0lMp(q, k) and M D (q' , k') = M D (q' , fc') ! ^= • (2-10) 

Here ul{p) is the left-handed neutrino spinor, 3 p = q — k, p' = q' — k', and po, p' are the 
time components of the corresponding 4-momenta. 

The quantity Sp a {x 2 , x\) in the second line of eq. (2.9) is the coordinate-space neutrino 
propagator in matter in the flavour basis. It is a matrix in both flavour space and spinor 
space, whereas the quantities Mp(q,k) and Mpt{q\k') are Dirac spinors. To simplify the 
notation, we have suppressed the corresponding spinor indices. 

The neutrino propagator in matter Sp a (x,x') satisfies the Schwinger-Dyson equation 

[ilpd* - MP R - M^P L ]S(x,x') - J d A x"Z(x,x")PS(x",x') = 5 A (x -x')-t, (2.11) 

where M is the neutrino mass matrix in the flavour basis, Ti(x,x') is the matter-induced 
neutrino self-energy, 1 is the unit matrix in the flavour space, and flavour indices are 
suppressed for simplicity. Eq. (2.11) (as well as eq. (2.13) below) applies to both Dirac 
and Major ana neutrino cases, provided that for Major ana neutrinos one uses the Feynman 
rules with propagators and vertices not containing explicitly the charge-conjugation matrix 
(see, e.g., [11, 12]). The operator P in (2.11) is defined as P = Pj_, for Dirac neutrinos 
and P = —75 for Majorana neutrinos. We will discuss the choice of the operator P in the 
Majorana neutrino case in Appendix A. 

For Dirac neutrinos, the mass matrix M is a general Nf x Nf matrix where Nf is 
the number of light neutrino species. Note that for M 7^ the presence of the term 
MPp + Pi rather than the usual mass term M in (2.11) is required by the hermiticity 
of the Lagrangian. For Majorana neutrinos, M = M T . 

3 We do not write flavour indices for the neutrino spinors because in the limit of ultra-relativistic neutrinos 
that we consider the spinors ul{p) correspond to essentially massless neutrinos. 
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The matter-induced neutrino self-energy £(x,x') is due to neutrino interaction with 
the particles of the medium through the exchange of and Z° bosons. For low energies 
of neutrinos and background particles, one can expand the propagators of the and 
Z° bosons in the inverse powers of their squared mass; the leading (neutrino energy and 
momentum independent) terms in these expansions yield 

E(x,x') ~Y, (x)8 4 (x- x') , where E (x) = ^V(x) , (2.12) 

and V fl (x) can be considered as an effective neutrino potential. Next (finite-order) terms 
in the expansions in 1/m^/ and ]~/m 2 z bring in some neutrino momentum dependence, 
which in the coordinate representation would result in the appearance of derivative terms 
in eq. (2.12); however, upon integration over x" in eq. (2.11) these terms would still lead 
to local terms in the self-energy S. Thus, in the case when only a finite number of terms 
in the expansion of the propagators of the intermediate bosons in inverse powers of their 
squared mass is kept, neutrino interaction with matter can be described by a local effective 
potential V{x). A Eq. (2.11) can then be rewritten as 

[l^ - V»{x)P)5 M - M M P R - M^P L ] S ja (x, x') = 5\x - x')5 Q „ , (2.13) 

where we have reinstated the flavour indices, whereas spinor indices are suppressed as 
before. For a non-relativistic medium of unpolarized particles, only the time component of 
V^{x) is essentially non-zero: V^{x) ~ V(x)5 fl Q. The solution of eq. (2.13) for the neutrino 
propagator in the case of Dirac neutrinos was given in [9]. It is reviewed in Appendix A, 
where also the Majorana neutrino case is considered. 

Let us now discuss the spinor structure of the neutrino propagator and of the plane- 
wave production and detection amplitudes. The propagator can be represented as 

S = ( S q LL S q LR ) , (2.14) 

where Sll = PlSPl, Slr = PlSPr, etc., are 2x2 block matrices, and we have now 
omitted the flavour indices. Since only left-handed neutrinos participate in weak interac- 
tions, we are interested only in the Li?-component of the propagator: —i{TisL(x)vL( x ')) = 
PlS(x, x')Pr = Slr(x,x'). It will be convenient for us to work in the chiral (i.e. Weyl) 
representation of the Dirac 7-matrices, in which 75 is diagonal. It can then be shown 
that for ultra-relativistic neutrinos with spin down along the 3rd spatial axis only the 
22-component of the 2x2 matrix Slr(x,x') is non-zero [9]. Likewise, from the Dirac 
equation it follows that the left-handed 2-component neutrino spinors in the momentum 
space are ul(p) = (0, v / 2po) T (see, e.g., [13], eq. (3.53)). From eq. (2.10) we then find 

Mp(q,k)=M P2 (q,k) and M D (q' , k') = M D2 (q' , k') , (2.15) 



4 Note that the term 'potential' is not very precise. Strictly speaking it applies solely to the case when 
only leading terms of expansions in powers of l/m%- and l/m| are retained. Otherwise, V(x) would 
depend on neutrino energy. This will pose no problem if V(x) is time-independent (i.e. V(x) = V(x)), so 
that different neutrino energy modes can be studied separately. It would be more correct to call V(x) local 
matter-induced neutrino self energy. We use the term 'effective potential' for brevity. 
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where the index 2 stands for the second (lower) components of the left-handed spinors 
Mp(q, k) and Mp,{q' , k'). On the other hand, we have 

MdSlrMd = M D2 {S LR ) 22 M P2 , (2.16) 

where we have taken into account that only the 22-component of Slr is different from zero. 
Denoting this component as S, from (2.15) and (2.16) we find 

M D S LR M D = M D S M P . (2.17) 

We can now rewrite eq. (2.9) as 

iA$;-( q ,k;q',k') = J d 4 Xl J d 4 x 2 M D (q',k')e-^'- k '^-^ 

x S Pl (x 2 , Xi)M P (q, k) e -i(Q-k)(xi-x P ) _ (2.18) 

where the integrand does not contain any quantities with spinor indices. 

It will be convenient for us to express the coordinate-space neutrino propagator Sp a {x 2 , x±) 
as a Fourier transform of the momentum-space one. We will be assuming that the matter- 
induced potential of neutrinos depends on the spatial coordinate x but is time-independent: 

= y M (x). 5 In this case the system under consideration possesses translation invariance 
in time but not in space; as a result, the coordinate-space propagator depends on the times 
ti and t 2 only through their difference, but on the spatial coordinates xi and x 2 separately: 
S(x 2 ,xi) = S(t 2 — ti;x2,xi). This constitutes an important difference as compared to the 
case of neutrino propagation in vacuum, where the coordinate-space neutrino propagator 
depends only on x 2 — x\ = (t 2 — t±; x 2 — Xx). The momentum-space neutrino propagator 
in non-uniform but static matter will therefore depend on one energy variable po = p' = E 
and two momenta: S = S(E; p', p). The coordinate-space neutrino propagator is related to 
the momentum-space one through the Fourier transformations with respect to the energy 
and both momenta: 

V( i2 - il ;x 2 ,x 1 ) = / f ^|^V(i?;p',p)e-^-^e-'x 2 - J px 1 . (2 . 19) 

Substituting this into eq. (2.9), going to the shifted integration variables x[ = x\ — xp and 
x 2 = x 2 — xp> and then using the obtained result in eq. (2.8), we find 

Apa = J g ^ s ( f^ 3 * D (E, p') ~S Pa {E- p', p) * P {E, p) e -*gfa-*p) e V*»-*P*r . 

(2.20) 

Here the functions <&p(E, p) and &d(E,p') are given by 

* P {E, p)=j d^x'^ J [dq] J [dk] / Pi (q, Q) f* Pf (h, K) e^- k ^M P (q, k) , (2.21) 
* D (E,p') = J d 4 x' 2 e-^'^ J[dq'} J[dk'}f D M,Q')fh f ^,K')e-W- k '>*M D {q',k'), 



5 This is a good approximation provided that the potential is nearly static over the time intervals of order 
At = a xv /v v , where a xv is the length of the neutrino wave packet and v v ~ 1 is the neutrino velocity. 
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where the 4- vectors p and p' are defined as p = (E, p), p' = (E, p'). 

The quantities §p(E,p) and <&£>(.E, p') are the amplitudes of the neutrino production 
and detection processes in which the external particles are described by the state vectors 
(2.5) and (2.6), while the produced and detected neutrino states are described by plane 
waves of 4- momenta p and p' , respectively. They are thus the probability amplitudes 
that the emitted and detected neutrinos have the corresponding 4-momenta, i.e. are the 
amplitudes of the momentum distribution functions of these neutrinos. In the limit of plane- 
wave external particles, &p(E, p) and &p>(E, p') are proportional to J- functions expressing 
the momentum conservation at neutrino production and detection. In the realistic case 
when the external particles are described by wave packets, the functions &p(E, p) and 
$£>(£', p') represent approximate conservation of mean momenta in the neutrino production 
and detection processes and are characterized by peaks of finite widths, with peak momenta 
being, respectively, P Q — K and P' « Q' — K' [10]. The width a p p of the peak of 
the function &p(E,p) depends on the momentum uncertainties of the particles taking 
part in neutrino production. It is dominated by the largest of these uncertainties: a p p ~ 
max{<7pj, &pf}. Quite analogously, the width u p d of the peak of the function &d(E, p) 
satisfies a v p> ~ max{o"p,j, 0\d/}- 

By the Heisenberg uncertainty relations, the momentum uncertainties at neutrino pro- 
duction and detection are related to the spatial localizations of the neutrino production 
and detection processes, a x p and a x D'- 

11 

c P p ~ , <?pD ~ • (2.22) 

&xP &xD 

2.2 Transition amplitude: a simplification 

Let us now proceed with the calculation of the transition amplitude. In the case of neutrino 
oscillations in vacuum, there exists a closed- form expression for the neutrino propagator, 
which in the momentum space depends on just one 4- momentum p. For the coordinate- 
space propagator one can then use an asymptotic expression at large baselines L given by 
the so-called Grimus-Stockinger theorem [14]. This leads to a considerable simplification of 
the expression for the amplitude of the overall neutrino production-propagation-detection 
process [10]. Unfortunately, for neutrino oscillations in matter with an arbitrary den- 
sity profile no closed- form expression for the neutrino propagator exists, and the Grimus- 
Stockinger theorem cannot be utilized. One therefore has to find another way to proceed 
with the computation. We shall show now that the calculations can be greatly simplified by 
making use of the fact that the momentum dependence of different factors in the integrand 
in eq. (2.20) has different character. 

Let us first note that the phase p'xp, — pxp of the momentum-dependent complex 
phase factor in (2.20) can be written as ^(p' — p)(x£> + xp) + |(p' + p)(xp, — xp). The 
first term here can be eliminated by the proper choice of the origin of the coordinate 
frame. The second term implies that the exponential factor in the integrand of (2.20) varies 
significantly when the momenta p, p' vary by |Ap|, |Ap'| ~ L , where L = |x£> — xp| is 
the baseline. Since L is a macroscopic distance, the phase factor in the integrand of (2.20) 
is a fast oscillating function of the momenta p and p'. At the same time, the neutrino 
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production and detection amplitudes <frp(E, p) and &d(E, p') are slowly varying functions 
of the momenta. Indeed, they change significantly when the corresponding momenta vary 
by |Ap| ~ a p p ~ l/VxP and |Ap'| ~ a p p> ~ l/cr x D- Because the sizes of localization 
regions of the neutrino production and detection regions are by far much smaller than the 
oscillation baselines of interest, 6 the amplitudes Qp(E, p) and $d(E, p') change very little 
over the momentum intervals over which the phase factor varies significantly. Therefore, 
these amplitudes can be pulled out of the momentum integrals at the values of momenta 
p = p* and p' = p*, where p* and p* are the central momenta of the regions which give 
the main contributions to the integrals over p and p', respectively. As a result, eq. (2.20) 
becomes 

Ap a = J ^ D (E,ti)$ P (E,p*)e- iE ^-^ J ( ^3 0js SpaiE; P', P) e**°-^ . 

(2.23) 

The last integral here is nothing but the neutrino propagator in the mixed energy-coordinate 
representation: 

% a (E;x Dl x P ) = J dre^ Q (r;x D ,x P ) = J ^ ^ Sp a (E; p', p) e ^^~^ , 

(2.24) 

where the last equality follows from eq. (2.19). Thus, we finally obtain 

/dE 
^d(E, p'J ®p(E, p*) e - iE ^~ tp ) Sp a (E; x D , x P ) . (2.25) 

Let us now discuss the propagator Sp a (E;x,x.'). It has been shown in [9] (see also 
Appendix A) that this quantity can be represented as 

i|i?||x— x'l 

Sp a (E; x, x') = -2E^- A{E- x, x') , (2.26) 

47T|X — X | 

where E > for neutrinos, E < for antineutrinos, and Fp a (E;x,x!) satisfies the equation 



—F 

dx 



w +y(x) 



F , where -^-=r-V. (2.27) 
dx 



Here r is the unit vector in the direction of the neutrino propagation: r = (x — x') /|x — x'| , 
and -4- as the directional derivative along f . The effective potential V(x) is related to the 
components of the neutrino potential in matter V^(x) through 

V = V° -v u -\~ ~V° -V 3 , (2.28) 

where v u is the neutrino velocity vector and V 3 is the component of V in the direction of 
neutrino propagation. The potential for antineutrinos is obtained from that for neutrinos by 
flipping the sign of the latter (except in CP-symmetric or nearly CP-symmetric media, see 



6 For instance, if a x p and a X D are of the order of interatomic distances and L ~ 1 km, then 
{a*P,(T*D}/L~ 1(T 13 . 
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Appendix B). Taking into account that MM^ = UM^W (which is valid in both the Dirac 
and Majorana neutrino cases), and that the neutrino potentials that enter in eqs. (2.27) 
and (1.3) are defined in the same way, we find that these equations coincide. Thus, the 
quantity Fp a satisfies the same equation as the amplitude of v a — > up oscillations in the 
standard approach to neutrino oscillations in matter. We summarize the expressions for 
the potential V{x) for neutrino propagation in various media in Appendix B. 

Let us now return to eq. (2.25). We have defined p* and p'^, as the momenta, small 
neighbourhoods of which give the main contributions to the integrals over p and p' in (2.24). 
How can one find these momenta? In the case of vacuum neutrino oscillations, the Grimus- 
Stockinger theorem tells us that due to a fast oscillating phase factor in the integrand of 
the Fourier-integral representation of the coordinate-space neutrino propagator, the neu- 
trino is forced to be on its mass shell. Hence, &p(E, p) and §p>(E,p') should also be the 
production and detection amplitudes for on-shell neutrinos. This is a simple consequence 
of the fact that particles propagating macroscopic distances are essentially on their mass 
shells. One can therefore expect that for neutrino oscillations in non-uniform matter, in 
the case when neutrinos propagate macroscopic distances the production and detection 
amplitudes &p(E, p) and &p)(E,p') should also be taken on the "in-matter mass shells" 
corresponding to the neutrino production and detection points, respectively. Here by the 
"in-matter mass shell" we mean that the neutrino energy and momentum at a fixed point 
with coordinate x should satisfy a dispersion relation that follows from the neutrino evolu- 
tion equation (2.27) with the effective potential V(x). Thus, we expect that the momenta 
p* and p'^ should satisfy the in-matter dispersion relations with the potentials V(xp) and 
V(xp>), respectively. A direct proof of this statement will be given in Appendix C. 

Let us now show that the in-matter dispersion relations are well defined for neutrino 
production and detection. Indeed, the sizes of the localization regions of the neutrino pro- 
duction and detection processes, a x p and ct x d, are small in comparison with the typical 
distances over which the neutrino potential in matter V(x) varies significantly. There- 
fore, to a very good accuracy one can consider the production and detection processes 
as occurring at constant densities given by the matter densities at, respectively, neutrino 
production and detection points xp and xj> 

Note that the transverse components of the neutrino momentum are extremely small, 
Ip_l|/.P ~ maxl^p,^}/! < 10~ 13 (see footnote 6), and so they can be safely neglected. 
As shown in Appendix C, the longitudinal components of the characteristic momenta p* 
and p^ satisfy eqs. (C18) and (C19). Since these are matrix equations, it is convenient to 
go to the basis where the matrix H is diagonal. For this purpose, we introduce the neutrino 
mixing matrix in matter according to 




(2.29) 



K 
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The states |i/#(x)) are thus the local eigenstates of H(x), which are called local matter 
eigenstates. Eq. (2.29) relates the neutrino flavour eigenstate basis to the basis of the 
local matter eigenstates |i/#-(x)), just like eq. (1.2) relates it to the mass eigenstate basis. 
In the limit of vanishing matter density the mixing matrix in matter U(x) goes to the 
vacuum mixing matrix U and the matter eigenstates go to the mass eigenstates. Note that 
eq. (2.29) merely describes a basis transformation in eq. (1.3); it does not necessarily define 
the matter-eigenstate content of the produced and detected neutrino flavour states. It only 
does so when the neutrino production and detection coherence conditions are satisfied. We 
will discuss this point in more detail in section 4. 

We also introduce the neutrino propagator in the matter eigenstate basis Sk'k{E\ x', x). 
According to eq. (2.29), it is related to the flavour-basis propagator Sp a {E\ x', x) through 

Sp a (E\ x , x) 

K,K< 

From (2.26) it follows that there is a similar relation between Fg a and the corresponding 
matter-eigenstate quantity Tk'K'- 

i^ a (£;x',x) = Upk'{^)Kk(^K'k(E;x! ,x) = [U{x')F{E;x! ,yi)U^{x)]p a . (2.32) 

K,K' 

Going in eqs. (C18) and (C19) to the matter-eigenstate basis allows one to immediately 
solve them with respect to the momenta p* and p*, which are the longitudinal components 
of p* and p' + . The results are given in eqs. (C20) and (C21). Since the transverse compo- 
nents of p* and p* essentially vanish, we conclude that the in-matter neutrino dispersion 
relations fully define the momenta that give main contributions to the integrals over p and 
p' in (2.24). Indeed, eqs. (C20)-(C22) imply p* = p^ and p* = p' K ,, where px and p' K , 
are the momenta of neutrino matter eigenstates at the production and detection points 
respectively. With this identification of p* and p*, eq. (2.25) can be rewritten as 




(2.33) 



This is the expression that we will be using in the following. 

3 Total rate of the process and the oscillation probability 

In the previous section we calculated the amplitude of the overall neutrino production, 
propagation and detection process. Our next goal is to calculate the probability of this 
process and then extract from it the oscillation probability. 

Let us first recall how the oscillation probability is determined from experimental 
data. Assume that, in an experiment, neutrinos of flavour a are emitted by a source, with 
the neutrino production rate and energy spectrum being T^ od and dT^ od {E)/dE. Let a 
detector sensitive to vp be situated at a distance L from the source, and the detection cross 
section be ap{E). The rate of the detection process is then 

rg = J dEjp(E)ap(E). (3.1) 
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Here jp{E) is the flux of vg at the detector site, which is given by 

, , i drl TOd (E) . . 

jp{E) = 4vtL2 dE P M E >*d,xp), (3.2) 

where P a p(E, x,d, xp) is the neutrino oscillation probability, and we assumed for simplicity 
that neutrino emission is isotropic. Substituting (3.2) into (3.1) yields the rate of the 
overall production-propagation-detection process: 

L a/3 = dE — — — - I dh — P a p(E,x D ,x P )crp(E) . (6.6) 

If the spectral density of the overall process rate dT t °^j(E)/dE is experimentally measured, 
one can find the oscillation probability by dividing this spectral density by the production 
rate, detection cross section and the geometric factor l/4irL 2 : 



dT™(E)/dE 
^[dTr d (E)/dE]^(E) 



Pap(E,x D ,x P ) = - _ pr :" — • (3.4) 



Notice that an important ingredient of this argument is the assumption that, for a fixed 
neutrino energy, the overall rate of the process factorizes into the production rate, os- 
cillation probability and detection cross section. If such a factorization turns out to be 
impossible, the very notion of the oscillation probability loses its sense, and one has to deal 
instead with the probability of the overall process. 

Now, we shall calculate the rate of the overall process in our QFT-based approach and 
try to present it in a form similar to (3.3), which would allow us to find the oscillation prob- 
ability. In doing so, we shall be closely following the treatment of the vacuum oscillations 
case in section 5.2 of ref. [10], to which we refer the reader for details. 

Let us first calculate the spectral density of the production rate dT^ od (E) / dE and the 
detection cross section ap(E). To simplify the calculation, we will be assuming that the 
neutrino emission and absorption processes are isotropic (relaxing this assumption would 
complicate the analysis but would not change the final result for the probability of neutrino 
oscillations) . This means that we can average the production and detection amplitudes over 
the direction of the incoming particles Pi and -D,, which amounts to averaging over the 
directions ofL = xd - xp. One can therefore define 3>p(-E, j>k) — J* ^gf-^p{E, psr), 
§ D (E,p' K ,) = J ^$ D (E,p' K ,). Applying the standard QFT rules, one then finds for the 
neutrino production and detection probabilities 

d 3 p 



p r * = J2\U aK \ 2 J ^S|*p(^P^)| 2 = EI^| 2 ^2 / dE\$ P (E, PK )\ 2 Ep K , 

(3.5) 

P$ e \E)=J2\UpK'\ 2 \*D(E,p' K ,)\ 2 -^, (3.6) 

K> N 

where Vjv is the normalization volume, and px, p^-/ are the energy-dependent momenta 
of neutrino matter eigenstates for V(x) = V(xp) and V(x) = V(x£>), respectively, which 
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are given by eqs. (C20) and (C21). In eqs. (3.5), (3.6) and in the following we use the 
shorthand notation 

UaK = UaK(*p) , U/3K' = U/3K'(xd) , (3.7) 



i.e. a and K,M,... will always refer to, respectively, the flavour index of the produced 
neutrino state and the indices of its matter-eigenstate components, whereas (3 and K', M' , ... 
will similarly refer to the flavour and mass-eigenstate components of the detected state. 

From eqs. (3.5) and (3.6) one can find the spectral density of the produced neutrino 
flux and the detection cross section [10]: 



dE T 2vH 



(3.8) 
(3.9) 



K' 



Here Np/Tq and Nd/Tq are flux-dependent normalization constants [10], which will drop 
out of the final result for the oscillation probability. 

Next, we need the rate of the overall neutrino production-propagation-detection pro- 
cess, which can be found by integrating the squared modulus of the amplitude of the process 
over the production and detection times tp and tp>. The time integrals can be reduced to 
the integrals over (tp + tp>) jl and T = tp> — tp. The first integration is trivial, whereas the 
second one leads to 

r$ = ^£ J dT|^ Q (T,x D ,x P )| 2 . (3.10) 
Substituting here the expression for the amplitude A/3 a (T, xp, xp) from (2.33), we find 

rffi = Np £ D * 2 / ^(m 2 E K K u m ,u aM ui M ^ D {E^' KI )^p{E, VK ) 

^ ' J K,K',M,M' 

j Pm)J^k'k(E', xp,, ~x.p)J-^i m (E\ xp,, xp) . 

(3.11) 



The quantity J : k'k(E;xd,xp) introduced here is related to Sk'k(E;xd,Xp) in the same 
way as Fk'k(E; Xp,, xp) is related to Sk'k(E; xp, xp) (see eq. (2.26)). The spectral density 
dT t °p(E)/dE is obtained from the right hand side of eq. (3.11) by removing the integration 
over E. 

By comparing eq. (3.11) with eqs. (3.8) and (3.9) it can be seen that the factorization of 
the rate of the overall process into the production rate, propagation (oscillation) probability 
and detection cross section as in eq. (3.3) is only possible if the production and detection 
amplitudes ^p, <l?p can be pulled out of the sum in (3.11). This, in turn, is allowed only 
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if the corresponding momenta of the matter eigenstates satisfy 

|PK - Pm\ < crpP , (3.12) 
\p'k> ~ Pm'I < a pD ■ (3.13) 

Indeed, under these conditions the factors &p(E,pk) and &d(E, p' K ,) are essentially in- 
dependent of the indices K and K'; one can therefore replace them, respectively, by the 
quantities <&p(E, p) and &p>(E, p') calculated at the mean momenta p and p' and pull them 
out of the sum. From eq. (3.11) we then find 



N P N D 1 2 /N|2 

-l$p(^,p)r i$d(^,p /ni - 



d£ T 2 (4vr) 2 L 2 1 

X X] ^aKUf3K'U a MUpM'^ r K'K{E]^D,^p)^ : M'Al(E;^D,^p) ■ 
K,K',M,M' 

(3.14) 

Likewise, under conditions (3.12) and (3.13) one can replace <l?p and as well as the 
factors px and 1/pk' m eqs. (3.8) and (3.9) by the corresponding quantities taken at the 
average values of the relevant momenta. They can then be pulled out of the sums, which 
yields 

ap(E) = ^\<Z> D (E,pi)\^. (3.16) 

Here we have used unitarity of the leptonic mixing matrix in matter U. Substituting these 
expressions, together with dT t °^(E)/dE from eq. (3.14), into (3.4), we arrive at 

Pa/3(E, x D,*p) = ^2 U* K tjfSK'U a MU^ M ,J : K'K(E]X D ,Xp)Tl I/M (E-,X D ,Xp) 
K,K' S M,M' 

= \Fp a (E;-x D ,xp)\ 2 , (3.17) 



where the flavour-basis function E/j a (E;xp>,xp) obeys eq. (2.27) with the boundary con- 
dition (A10). Here the factors \$>p(E, p)| 2 |$£>(i£, p')| 2 in the numerator and denominator 
have canceled out, leaving us with the oscillation probability that is independent of the 
neutrino production and detection processes. In deriving eq. (3.17) we have also canceled 
p and p /_1 in the product (dr a rod (E) / 'dE) x crp{E) in the denominator. This is justified 
because the mean neutrino momenta at production and detection coincide to a very good 
accuracy under the conditions Am 2 / (2E) <C E, \V\ <C \E\, which we assume to be satisfied 
throughout this paper. 



7 While conditions (3.12) and (3.13) ensure the production and detection coherence, they say nothing 
about another possible source of decoherence - separation of neutrino wave packets at long enough distances 
L > Lcoh due to the difference of the group velocities of different neutrino mass eigenstates. This is related 
to the fact that a fixed neutrino energy corresponds to the stationary situation, when the coherence length 
icoh — > oo. The finite coherence length is recovered upon the integration over energy in eq. (3.11) [15]. 
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Thus, we have found that under conditions (3.12) and (3.13) the oscillation prob- 
ability can be sensibly defined and can be extracted from the rate of the overall neu- 
trino production-propagation-detection process. Since the matrix F is unitary, 8 the re- 
sulting oscillation probability (3.17) obeys the unitarity conditions P a p{E, X£>, xp) = 
X^a Pap{E, X£>, xp) = 1, i.e. is properly normalized. If conditions (3.12) and (3.13) are 
not fulfilled, the oscillation probability cannot be defined, and flavour transitions should 
instead be described by the rate of the overall neutrino production-propagation-detection 
process (3.11). 

Eqs. (3.12) and (3.13) are actually the conditions of coherent neutrino production 
and detection: their fulfilment ensures that the production and detection processes cannot 
distinguish between different neutrino matter eigenstates, so that these eigenstates are pro- 
duced and detected coherently. If these conditions are violated, i.e. if either |p/<— Pa/| ^ CpP 
or \p' K , — p' M ,| > a p D, the differences of momenta of different matter eigenstates will exceed 
the momentum widths of the corresponding momentum distribution amplitudes, <I>p or <I>p>. 
In that case the overlap of the amplitudes corresponding to different matter eigenstates 
will be suppressed, leading to a quenching of the interference terms in expression (3.11) for 
the probability of the overall process. Note that the momentum uncertainties due to the 
localization of the neutrino production and detection processes, a p p and cr p D, are usually 
much smaller than the neutrino momentum itself; therefore, conditions in eqs. (3.12) and 
(3.13) are much stronger than the conditions \pk~ Pm\ "C Pk,Pm, \p'k'~ Pm'I ^ Pp>Pm'i 
which follow automatically from Am 2 /(2E) <C E, \V\ <C \E\. 

4 The amplitude of the overall process, the oscillation amplitude and 
their evolution equations 

We have demonstrated in the previous section that in the case when neutrinos are ultra- 
relativistic, the matter-induced neutrino potential satisfies |^(x)| <C \E\, and in addition 
the conditions of coherent neutrino production and detection (3.12) and (3.13) are fulfilled, 
the oscillation probability can be sensibly defined and can be extracted from the rate of the 
overall neutrino production-propagation-detection process. The resulting expression for the 
oscillation probability in eq. (3.17) is simply given by the squared modulus of Fp a , which 
therefore can be interpreted as the oscillation amplitude. As we have already discussed, Fa a 
satisfies the evolution equation (2.27) (which coincides with eq. (1.3)), supplemented by the 
boundary condition (A10). Thus, in the case when the coherence conditions for neutrino 
production and detection are satisfied, the standard approach to neutrino oscillations in 
matter based on the Schrodinger-like evolution equation (1.3) is justified. 

Let us now discuss the amplitude of the overall neutrino production-propagation- 
detection process. Does it satisfy an evolution equation similar to (1-3)? Consider first 
the case of vacuum neutrino oscillations. The neutrino production and detection coherence 
conditions now read 

VpP,crpD- (4.1) 

s This follows from the fact that F satisfies the Schrodinger-like equation (2.27) with the Hermitian 
effective Hamiltonian, supplemented the boundary condition (A10). 
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Here pj = (E 2 — m 2 ) 1 / 2 ~ E — m 2 /(2E) is the momentum of the jth neutrino mass 
eigenstate of energy E. If these conditions are satisfied, the oscillation amplitude can be 
defined, and it coincides with the standard amplitude of neutrino oscillations in vacuum: 



Am . , 

jk r 



[A~(E, x)]p a = £ . U* aj U pj e-^ x . (4.2) 



The amplitude (4.2) satisfies the Schrodinger-like equation 



d 



dx 



.Am 



2 



2E 



i-A°^{E,x) = U—U^ A°T c (E,x) (4.3) 



with the boundary condition [A° s ^ c (E,0)]p a = 5p a . 

Let us now examine the probability of the overall neutrino production-propagation- 
detection process in vacuum, without assuming anything about coherence of neutrino pro- 
duction and detection. This probability can be written as [10] 

T%{x) = ^ J ^(2E) 2 \A^ c (E,x), a \ 2 , (4.4) 

where the quantity 

A^(E,x)p a = Y,KF^p(E,Vi)®D{E,Viy^> = {U[$ P <!> D e iA ntf}p a (4.5) 
j 

can be considered as the amplitude of the overall process. It is eqs. (4.4) and (4.5) that 
have to be used to describe neutrino flavour transitions in vacuum in the case when the 
coherence condition in eq. (4.1) are violated. From eq. (4.5) it is easy to find that the 
amplitude of the overall process in vacuum A^.(E,x) satisfies the same evolution equa- 
tion as the oscillation amplitude. Indeed, differentiating (4.5) we obtain i{d/dx)A t °^ c = 
U[<5> P $ D e iApx (-Ap)]U^ = U(-Ap)UiU[& P $ D J Apx ]rf = U(-Ap)rfA™, which coin- 
cides with (4.3). Crucial to this derivation was the point that all the factors in the square 
brackets are diagonal and therefore commute with each other. 

Although A° s £.(E,x) and A^ C (E, x) satisfy the same evolution equations, the bound- 
ary conditions that they obey are different. As was mentioned above, for the oscillation 
amplitude it is the standard condition [A^.(E, 0)]/3 a = Sp a ; at the same time, for the 
overall amplitude the boundary condition is [A^.(E, 0)]p a = {U$p$dU^} p a , as can be 
immediately seen from eq. (4.5). Obviously, the solution of one and the same eq. (4.3) with 
two different boundary conditions are different. 

Now let us return to neutrino oscillations in matter. The rate of the overall pro- 
cess (3.11) can be cast in the same form as in eq. (4.4), but with the vacuum amplitude 
A^.(E,x)p a replaced by 

A%(E,x,xo) = J2 Kk(^)Up KI {^p(E, Vk )^ d (E, v ' ki )Pk'k{E-^^) 



K,K' 



{U(x)$ D f<S> P W(x )}p a . (4.6) 
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This expression has a simple physical interpretation: the factor U^(xq) projects the ini- 
tial flavour-eigenstate neutrino u a onto the matter eigenstate basis, &p(E, px) are the 
amplitudes of production at the point xo of various matter eigenstates that compose v a , 
^(-E^XjXo) describes the propagation of these matter eigenstates to the point x (including 
the transitions between them), $£>(£", p' K i) are the detection amplitudes of neutrino matter 
eigenstates at the point x, and finally U(x) projects the amplitude back from the matter 
eigenstate basis to the flavour basis. 

Consider the case when both the coherence conditions (3.12) and (3.13) are violated, 
so that the amplitudes $p and $d cannot be pulled out of the sum in (4.6). Does the 
amplitude of the overall process A satisfy the same evolution equation as the quantity F, 
as it is the case for neutrino oscillations in vacuum? By differentiating eq. (4.6) with respect 
to x, 9 we immediately find that in general this is not the case. The reason for this is that, 
unlike in the case of vacuum oscillations, the matrix T{E\ x, xo) is not diagonal. Actually, 
for neutrinos moving in non-uniform matter the neutrino propagator is not diagonal in any 
basis. This comes about because the effective Hamiltonian H(x) cannot be diagonalized 
by one and the same unitary transformation for all values of x. The only exception is the 
special case of adiabatic neutrino evolution, when the propagator is diagonal in the matter 
eigenstate basis. In this case, by differentiating (4.6) with respect to x one can make sure 
that the oscillation amplitude satisfies the standard evolution equation (1.3) (though with 
a non-standard boundary condition). The proof is very similar to the one in the case of 
vacuum neutrino oscillations and is given in Appendix D. 

What happens in the situations when one of the coherence condition (3.12), (3.13) is 
satisfied, while the other is not? To answer this question, it will be convenient for us to 
rewrite eq. (4.6) in the form 

Af*(E,-x,xo) = {c7(x)cI> D L>t(x)F(^;x,xo)c>(x )$p[7t(xo)} /3Q , (4.7) 

where we have used eq. (2.32). This expression admits a simple interpretation similar to 
that of eq. (4.6) (see below). 

Consider first the case when the detection coherence condition (3.13) is satisfied, but 
the production coherence condition (3.12) is violated. In this case one can replace the 
factors $>d(E,p' k ,) in eq. (4.7) by the one taken at the mean momentum p', which yields 

A^(E,x,xo) = <S> D (E,pi){F(E;x,x )U(x )$pU\x )}p a - (4.8) 

From the fact that only the first factor in the curly brackets here depends on x, it imme- 
diately follows that in this case the amplitude A tot (E, x, xo) satisfies the same equation 
as F(E;x,xq) does, i.e. eq. (2.27). The boundary condition for the overall amplitude is, 
however, different: from eq. (4.8) we find 

Ag(S,x,xo)| x -.xo = ^D(^,p , ){^(xo)^P^ t (x )} /3a . (4.9) 

9 Recall that d/dx here is understood not as the derivative with respect to |x|, but as a directional 
derivative along r = x — xo, see eq. (2.27). 
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Now let us consider the opposite case when the production coherence condition (3.12) 
is satisfied, but the detection coherence condition (3.13) is not. Then from (4.7) we find 

A^(E, x, x ) = $ P (E, p){l) (x)$ d U\-k) F(E; x, x )}^ . (4.10) 

This expression contains, in addition to F(E;x, xo), two more x-dependent factors, U(x.) 
and C/(x)t; it can be readily seen that the amplitude .A tot (-E, x, xo) does not satisfy the 
same equation as F(E;x, xo) in this case. 

Thus, we found some disparity between the production and detection processes: if 
neutrino detection is coherent but the production process is incoherent, the amplitude 
of the overall process obeys the standard evolution equation (1.3), while in the opposite 
situation it does not. The reason for this asymmetry is that we assume the neutrino 
production coordinate to be fixed and consider the evolution of the amplitude with the 
coordinate of the neutrino detection point. If the detection process is incoherent, the 
flavour-eigenstate detection amplitude f/(x)<3?£)C7^(x) is coordinate-dependent, and the x- 
dependence of A tot (E, x, xo) is different from that of F(E;x,xq). Therefore the amplitude 
A tot (E, x, xo) does not satisfy eq. (1.3). 

How can one understand the above results in physical terms? Consider the matter 
eigenstate content of the initially produced neutrino state u a . The probability amplitude 
that the initial flavour state contains the matter eigenstate vk is given by &p(E, pk)U* k (x.q) 
It differs from the naively expected factor U* k (xq) that would follow from eq. (2.29) by the 
presence of the independent amplitude of vk production Qp(E, Pk)- In general, eq. (2.29) 
should actually be considered as the definition of the matter eigenstate basis rather than 
a relation giving the matter-eigenstate composition of the flavour neutrino state, which is 
process dependent. This comes about because eq. (2.29) describes the basis transforma- 
tion in the evolution equation (1.3) which ignores the coherence issues. If the production 
coherence condition (3.12) is satisfied, all the amplitudes <&p{E,px) can to a very good 
accuracy be replaced by a common factor §p(E, p). In this case the relative weights of 
different matter eigenstates in u a are given by | \J a K (xo | 2 , i-e. eq. (2.29) does give the 
matter-eigenstate content of v a . If, on the contrary, condition (3.12) is strongly violated, 
different matter eigenstates will be produced incoherently. Indeed, the squared modulus 
of the overall amplitude contains terms proportional to &p(E,pk)&*p(E,pm)] for K ^ M 
these are the interference terms. If \pm — Pk\ is large compared to the momentum width 
a p p of the amplitude <3>p, the quantities $>p(E,pk) and <3?p(£,pj\/) will have little overlap. 
In this case the interference terms are strongly suppressed, which means that vk and vm 
are emitted incoherently. 

The initially produced neutrino state can then be evolved from xo to x by F{E; x, xo), 
as in eq. (4.6). Alternatively, one can project the initial state onto the flavour basis and 
evolve it with F(E;x,xq), as in eq. (4.7). The evolved neutrino state is then absorbed in 
the detection process. Each of the matter eigenstate components fjc' of the evolved state 
is detected with its own amplitude &d(E, p' K i), so that in general the matter-eigenstate 
composition of the detected state is not given by eq. (2.29). Therefore the amplitude of the 
overall production-propagation-detection process does not satisfy the standard evolution 
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equation (1.3). However, if all &d(E, p' K ,) are to a good accuracy equal to each other (which 
is the case when the detection coherence condition (3.13) is fulfilled), the detection efficiency 
is essentially the same for all matter eigenstates, so that the detected flavour state is indeed 
related to the matter eigenstates by eq. (2.29). In this case the standard evolution equation 
(1.3) applies. As follows from the above discussion, this holds irrespectively of whether or 
not the production coherence condition (3.12) is obeyed. The latter just determines the 
initial state of neutrino evolution. 

In brief, if the matter-eigenstate composition of the evolving neutrino state is described 
by eq. (2.29), the amplitude of the overall process evolves according to eq. (1.3). Otherwise, 
eq. (1.3) does not apply, the only exception being the case of adiabatic neutrino evolution. 

There is an important remark that has to be added to the above discussion. We have 
found that in the case when neutrino detection is coherent while its production is not the 
amplitude of the overall process still satisfies evolution equation (1.3). However, even in 
this case the standard approach to neutrino oscillations in non-uniform matter has to be 
modified. This follows from the fact that the amplitude of the overall process does not 
factorize into the production, oscillation and detection amplitudes in this case; only the 
detection amplitude can be factored out. In such a situation one has to deal with the 
probability of the overall process, described by eq. (3.11). Alternatively, one can employ 
eq. (4.4) with the vacuum amplitude J\^ C (E, x)p a replaced by A^(E, x), where A^(E, x) 
satisfies eq. (1.3) with the boundary condition (4.9). 

5 Discussion and summary 

In this paper we have considered neutrino oscillations in non-uniform matter in the frame- 
work of QFT. We treated neutrino production, propagation and detection as a single pro- 
cess, described by the Feynman diagram of fig. 1, with neutrino in the intermediate state 
described by a propagator. We found that under certain conditions (which are satisfied in 
most cases of practical interest) the oscillation probability can be sensibly defined. 

We have demonstrated that when the conditions for the existence of the oscillation 
probability are fulfilled, this probability is given by eq. (3.17). The oscillation amplitude 
in this case coincides with the function Fp a (E; x, xo) that is simply related to the neutrino 
propagator in matter. This function satisfies the usual Schrodinger-like evolution equation 
(2.27), provided that matter density is sufficiently smooth, so that condition (A15) is 
satisfied. Hence, in this case the standard approach to neutrino oscillations in non-uniform 
matter is justified. We thus presented here a consistent derivation of the standard evolution 
equation, and found the conditions under which it is valid. Let us summarize here these 
condition once again: 

(i) Neutrinos are ultra-relativistic, so that 4rgr- <C E. 

(ii) The effective matter-induced potential of neutrinos depends on the coordinate but 
does not vary with time, i.e. V = V(x). 

(iii) The potential V(x) is small compared to the mean neutrino energy: |V(x)| <C E. 
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(iv) In the neutrino production and detection regions, matter density (and so the potential 
V(x)) is nearly constant. That is, V(x) varies little over the distances of order of the 
sizes of the spatial localization regions of neutrino production and detection, a x p and 
a x D- in other words, |V'/V| <C mm{a p p, ct p d}, where o p p ~ l/<r x p and ct p d ~ l/cr x D 
are the momentum uncertainties at neutrino production and detection, respectively. 

(v) Neutrino emission and detection are coherent, i.e. the conditions |p# — Pa/ I *C o p p 
and \p' K , — p' M ,\ <C (t p d are satisfied. 

In addition, when deriving eq. (2.27) we had to assume that the potential V(x) varies 
little over the distances of order of the neutrino de Broglie wavelength. However, since the 
momentum uncertainties at neutrino production and detection satisfy {cr p p, (T p d} *C p, this 
condition is superseded by the one in (iv), provided that the condition |V'(x)/V(x)| <C 
min{ (7 p p, a p r>} of point (iv) is fulfilled for all x along the neutrino trajectory and not only 
in the production and detection regions. 

Conditions (i)-(v) ensure that the oscillation probability can be sensibly defined and 
can be extracted from the probability of the overall neutrino production-propagation- 
detection process. Condition (iii) allows to simplify significantly the equation for F(E; x, Xq) 
and reduce it to the form (A20). Condition (ii) simplifies the consideration, but in fact 
is not necessary. It is enough to assume that the neutrino wave packets are sufficiently 
short, so that the potential is nearly constant in space and time over the distances of 
order a xu and times ~ cr xu /v u , where a xv is the length of the neutrino wave packet and 
v v is its group velocity. Note that this assumption is related to condition (iv) because 
°xv % max { (r ppi°'pi)}- Under the requirement that the potential V(x) vary very little over 
the distances ~ cr X y and times ~ ^xul^u the neutrino will "feel" a well defined potential 
along its path. If, in addition, o xv is small in comparison with the oscillation length and 
the baseline L, then one can consider neutrinos as pointlike particles. In this case, the 
potential can vary both in space and time, but at any point x on the neutrino trajectory 
only the value of the potential at the time t satisfying x = v u t will play a role, so that 

y(t,x) = y(|x|/K|,x) = y(x). 

As was discussed in section 3, the coherent neutrino production and detection condi- 
tions (3.12) and (3.13) are crucial for the possibility to define the oscillation probability as 
a production- and detection-independent quantity. If these conditions are not obeyed, one 
would have to deal instead with the rate of the overall neutrino production-propagation- 
detection process (3.11). The quantity J- that enters into this equation is related to F 
by eq. (2.32), while F should be found as the solution of eq. (2.27) with the boundary 
condition (A10). Flavour transitions are then, in general, not directly described by the 
standard neutrino evolution equation in matter (1.3). There are, however, exceptions from 
this rule. First, if the detection coherence condition (3.12) is satisfied, the amplitude of the 
overall neutrino production-propagation-detection process satisfies the standard evolution 
equation (1.3), supplemented by the boundary condition (4.9). This takes place even if the 
production coherence condition (3.12) is not obeyed and so the oscillation amplitude can- 
not be defined. Second, as shown in Appendix D, in the special case of adiabatic neutrino 
propagation the amplitude of the overall process satisfies the standard evolution equation 
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(1.3) even when both neutrino production and detection processes are not coherent. The 
boundary condition for the oscillation amplitude is given in this case by eq. (D5). 

Are there any situations in which the coherence conditions for neutrino production or 
detection (3.12), (3.13) are violated and therefore the oscillation amplitude satisfying the 
standard evolution equation (1.3) cannot be defined? As we shall see, this may only be 
possible for large values of the neutrino mass squared differences Am 2 , which would imply 
the existence of relatively heavy sterile neutrino states. 

Production and detection coherence conditions (3.12) and (3.13) actually require that 
the neutrino production and detection regions be small in comparison with the neutrino 
oscillation length (they can therefore be also called the localization conditions). Let us 
consider for simplicity a 2-flavour oscillation problem and discuss first neutrino production 
coherence. The production coherence condition (3.12) can then be written as 



' cos 26 -V(x P ) J +(^J sin 2 29 < a pP , (5.1) 



2E 



where 9q is the mixing angle in vacuum. Let us first consider the case when the neutrino 
potential at the production point dominates over the neutrino kinetic energy difference, 
i.e. |y(xp)| > Am 2 /{2E). Production coherence condition (5.1) is then violated when 

G F N e (x P ) > a pP > l/a xP . (5.2) 

Assume that the mean distance between the particles of the matter in the neutrino produc- 
tion region is r^. Then we have N e (x P ) ~ 1/Vjj, a xP ^ r 0i an d eq. (5.2) requires r 2 . < Gf, 
or ro < 6 • 10~ 17 cm. This corresponds to extremely high densities, exceeding the nuclear 
density by about ten orders of magnitude. Such densities are only attainable in the very 
early universe, when neutrino oscillations are irrelevant. 

Next, let us consider the opposite situation, |l/(xp)| <C Am 2 /{2E). Production coher- 
ence condition (5.1) is then violated if 

Am 2 

-fif ~ °vP ■ ( 5 - 3 ) 

Consider, e.g., neutrinos produced in an accelerator experiment in decays of pions of speed 
v n inside a decay tunnel of length l p . It has been shown in [33, 34] that in this case the 
production coherence condition is violated when 

Atti 2 Am 2 

—i p >i, (r/ P K«i) ; ^et v -> 1 (ry^»i), (5.4) 

where V is the pion decay width in the laboratory frame. In the case of relatively short 
decay tunnels (l p <C Idecay = v^/T) condition (5.4) yields l p > 2E/Am 2 = 1^/2tt, where 
/qsc is the vacuum oscillation length. Thus, in this case production coherence is violated 
when the length of the decay tunnel is comparable with the neutrino oscillation length. 

The opposite case of relatively long decay tunnels, l p 3> Idecay = v^/T, is, however, 
of more practical interest, since in this case most pions decay before being absorbed by 
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the wall at the end of the tunnel. In this case we have to use the second inequality in 
(5.4), which yields Am 2 > 1 eV 2 . Such values of Am 2 are currently widely discussed in 
connection with possible existence of light sterile neutrinos [35, 36]. 

Let us now briefly discuss possible detection coherence violation. As follows from 
our discussion above, one can concentrate on the case |V(x£»)| < Am 2 /(2E). Detection 
coherence condition (3.13) is then violated provided that 

Am 2 . . 1 . , 

^ttt £ VpD > , (5.5) 

IE a x o 

similarly to (5.3). Let the average distance between the particles in the detector be r$. 
Then <j x d < ro, and condition (5.5) requires Am 2 /(2E) > r Q 1 - For matter of normal 
density tq ~ 10~ 9 cm, and for neutrinos in the MeV range we find that condition (5.5) 
requires Am 2 > (100 keV) 2 . 

To summarize, we presented a consistent treatment of neutrino oscillations in non- 
inform matter within a QFT framework. We have found that the oscillation amplitude 
can be sensibly defined and can be extracted from the amplitude of the overall neutrino 
production-propagation-detection process if neutrinos are ultra-relativistic, matter density 
varies little over the distances of order of the sizes of the production and detection regions 
of individual neutrinos, and the neutrino production and detection processes are coher- 
ent. By the latter we mean that different matter eigenstates composing the flavour states 
are emitted and absorbed coherently. In this case the oscillation amplitude satisfies the 
standard evolution equation (1.3). Otherwise one has to consider instead the probability 
of the overall process, given in eq. (3.11). Production coherence can be violated e.g. in 
the case of neutrinos produced in decays of free pions provided that sterile neutrinos with 
Am 2 > 1 eV 2 exist and this mass squared difference plays a role in the flavour transitions 
of interest. For detection processes in matter of normal density (a few g/cm 3 ) one can 
expect coherence violation for Am 2 / (2E) > (100 keV) 2 /MeV. 

The authors are grateful to Georg Raffelt and Alexei Smirnov for very useful discus- 
sions. 



Appendix A: Neutrino propagator in non-uniform 
matter 

Here we briefly describe the calculation of the neutrino propagator in non-uniform matter 
in the Dirac and Major ana neutrino cases. In the Dirac case our treatment closely follows 
that of [9], the main difference being that we allow the neutrino mass matrix M to be an 
arbitrary non-singular matrix, whereas in [9] it was assumed to be hermitian. 

The coordinate-space neutrino propagator in matter satisfies eq. (2.13). We assume 
that the matter-induced neutrino potential is the function of the coordinate x along 
the neutrino trajectory but is time independent: V 1 * = V*(x). The neutrino propaga- 
tor Sp a (x,x') then depends on the times t and t' only through their difference, but on 
the spatial coordinates x and x' separately: Sp a (x,x') = Sp a (t — £';x,x'). It is conve- 
nient to introduce the neutrino propagator in the mixed energy-coordinate representation 
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Sp a (E;x,x'), which is related to Sp a (t — t';x,x') through 

Sp a (t-t';x,x') = J ^e- iE WSp a (E;x,x'). 



(Al) 



The inverse transformation is given by the first equality in eq. (2.24). 

From now on, we will distinguish between the Dirac and Majorana neutrino cases. 

A.l Dirac neutrino propagator 

In this case one has to set P = Pl in eq. (2.13). Omitting the flavour indices to simplify 
the notation and writing S(E;x,x') in the block-matrix form (2.14), from eq. (2.13) we 
find 



-Aft 

E - ia ■ V - V° - V • 



E + ia ■ V\ / S LL S L r 



-M 



Srl Srr 



<5 3 (x-x') 



1 
1 



(A2) 



We will only need the Sir block matrix of the neutrino propagator. From (A2) we obtain 
a system of two coupled equations for Slr and Srr: 

-M^S LR (E; x, x') + (E + ia • V)Srr(E; x, x') = , (A3) 
[(£ _ ia . V ) - y°(x) - V(x) • a}S LR (E;^') - MS RR (E;x,x') = 5 3 (x - x') . (A4) 

Next, define J(E;x,x') = (Aft)" 1 S RR (E;x,x'). Eq. (A3) then gives 

S LR (E; x, x') = (E + ia ■ V)J(E; x, x') . (A5) 

Substituting this into eq. (A4), we obtain the equation for J(E;x,x'): 

{E 2 + V 2 - MM 1 " - EV°(x) - iV(x) - V -ia - [F°(x)V - iEV(x) + iV(x) x V]} 

x J(E;x,x') = <5 3 (x-x'). 

(A6) 



We shall seek the solution of this equation in the form 

„i|E||x-x'| 



With this ansatz, 



VJ 



J(E;x,x') 



2 \E\e i \ E U x - x '\ 



47r|x — x 



-F(E;x,x') 



(V 2 + E 2 )J = 5 3 (x-x')F 



4tt\x — x'| 

2\E\e i ^ x -*'\ 
47r|x — x'l 



1 r „ 1 „ f 

— F + — VF - — -n 

2 2E 2\E x-x' 



(A7) 



(A8) 



2 LB 



■V 2 F + i(r ■ VF) 



\E\ x — x 



-(f-VF) 



(A9) 

where f = (x — x')/|x — x'|. Requiring the first term on the right hand side of (A9) to 
cancel the 5-function in eq. (A6) gives the boundary condition for F: 



i^ Q (f?;x,x') 



(AW) 
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where we have restored the flavour indices a and f3. Since in neutrino oscillation experi- 
ments one deals with macroscopic distances, we are primarily interested in well separated 
x and x'. This means \E\\x — x'\ 3> 1, so that the last terms in eqs. (A9) and (A8) can be 
neglected. Eq. (A6) then becomes 

i(r ■ VF) + 77^V 2 F - -L-[MAf1 + EV° - \E\(t • V) 
2 \h/\ ^\E\ 

- a ■ (V°\E\r -EV + i\E\V x t)]F + O M ^ — -j = , (All) 

Since in all situations of practical interest matter-induced neutrino potentials are by far 
much smaller than neutrino energy, 

\V\ < \E\ , (A12) 

the last term in (All) can be neglected in comparison with the first term. Choosing the 
z-axis of the coordinate system along f , one can then rewrite eq. (All) as 

i(f • VF) + V 2 F - ^D(E, x)F = , (A13) 

where 

' MMt + (E - \E\)(V° + V 3 ) (E + \EDiV 1 - iV 2 ) 
D{E,x) =1 I . (A14) 

{E -\E\){V l + iV 2 ) MM^ + (E + \E\)(V° -V 3 ) 

Let us now distinguish three cases: (1) |VF| 3> eF, where e is the largest eigenvalue of 
the matrix D(E,x)/2\E\; (2) |VF| < eF; and (3) \VF\ ~ eF. The first case (in which 
the third term in eq. (A13) can be neglected) is of no interest to us because it corresponds 
to the kinematic region in which neutrinos essentially do not oscillate. In the second 
case the three terms in eq. (A13) cannot balance each other, i.e. this equation cannot 
be satisfied. This immediately follows from (A12) and the condition Am 2 <C E 2 , where 
Am 2 = maxjAm-J. 10 Thus, the only case of interest to us is the third one. It is easy to 
see that in this case the second term in eq. (A13) is negligibly small compared to the other 
two and so can be omitted provided that neutrinos are relativistic, the components of the 
neutrino potential V^(x) satisfy eq. (A12), and in addition 

-^H«|£|- (A15) 

Note that this condition requires that the potential change little over the distances of order 
of the neutrino de Broglie wavelength. Under the above conditions eq. (A13) reduces to 

i(f-VF) r 1 — D(E, x)F = 0. (A16) 

2 \E\ 



10 



Obviously, only mass squared differences and not the absolute neutrino masses play a role in neutrino 
oscillations. Technically, this can be proven by subtracting from D(E,x.) the matrix mf [ ^ J, where mi 



1 

is any neutrino mass eigenvalue and 1 is the unit matrix in the flavour space, and rephasing F accordingly. 
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Let us now concentrate on the case of neutrinos, E > (the antineutrino case can be 
studied similarly). From |V-F| ~ eF we find that the second term on the right hand side 
of eq. (A8) is much smaller than the first one; we have already established that the third 
term in this equation is negligible. Thus, {E + icr ■ V) J ~ E(l — cr ■ r)J = E(l — u^)J . 
From eq. (A5) we then find 

(S LR ) 2 2 = 2EJ 22 , (A17) 

with all the other spinor components of Slr being zero. This result plays an important 
role in our calculations, since we need to deal only with one component of the neutrino 
propagator, and this simplifies our consideration significantly. Eq. (A7) relating Slr and 
F then implies that the only relevant spinor component of F is F^i- 
Next, we note that for E > eq. (A14) can be rewritten as 

/ MM* 2E(V 1 - iV 2 ) \ 
D{E,x) = . (A18) 

\ MM^ + 2E(V° - V 3 ) I 

The fact that D 2 ± = means, in particular, that the equation for the spinor component 
F 22 in (A16) decouples, i.e. does not contain any other components of F. 

Denoting the 22-components of Slr, J and F as (Slr)22 = S, J22 = J and F22 = F, 
we finally obtain from (A17), (A7), (A16) and (A18) 

e t|E||x-x'| . 

S LR (E; x, x') = -2E— -F(E; x, x') , (A19) 

47T|X — X I 

where F satisfies the Shrodinger-like equation 



dx 



MM] V, s 
+ F(x) 



F . (A20) 



2\E\ 

Here ^ as the directional derivative along f : 4- = r • V. The potential V(x) is defined as 

F(x) = T/°(x) - F 3 (x) ~ V°(x) - »„-V(x) , (A21) 

where v v {\v u \ ~ 1) is the neutrino velocity vector, and the last equality is valid for an 
arbitrary orientation of the coordinate axes. Eq. (A20) actually describes both the neutrino 
and antineutrino cases; the potential for antineutrinos is obtained from that for neutrinos 
by flipping the sign of the latter (except in CP-symmetric or nearly CP-symmetric media, 
see Appendix B). 11 

Note that the factor ^^j^L^ 1 in expression (A19) for Slr is a fast varying func- 
tion of the coordinates, which changes significantly over distances of order of the neu- 
trino de Broglie wavelength E , whereas eq. (A20) actually means that the factor F 



n The antineutrino case is studied quite analogously to the neutrino one. In that case one has to replace 
E — > — E, f — > — r. The only non-vanishing spinor components of Slr and J are S(lr)u and Jn, and the 
only relevant component of F is Fn, which satisfies the same eq. (A20) with the potential being negative of 
the neutrino potential (except in media with equal or almost equal numbers of particles and antiparticles) . 
The right-handed antineutrino spinors in the momentum space are vr(j>) = (y / 2po, 0) T , i.e. only their 
upper components are non-zero. 
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is a slowly varying function of x, which changes significantly over the distances of order 
min{£'/Am 2 , |U°| _1 , |V| -1 }, i.e. of order of neutrino oscillation length in matter. 

For a known matter-induced potential, eqs. (A19) and (A20) together with the bound- 
ary condition (A10) fully determine the neutrino propagator Slr(E,x,x'). 

A. 2 Majorana neutrino propagator 

Recall that for Majorana neutrinos we use the Feynman rules in which propagators and 
vertices do not contain explicitly the charge-conjugation matrix [11, 12]. Let is first discuss 
the choice P = —75 in eqs. (2.11) and (2.13). For Majorana neutrinos the 4-component 
field v = vl + ur can be written as v = ul + (vl) c , where the superscript c means 
charge conjugation. In other words, in this case right-handed neutrinos are antiparticles of 
left-handed ones, and so they participate in the standard weak interactions. The matter- 
induced potential enters the equations of motion of the right-handed and left-handed 
fields with opposite signs. The choice P = —75 in eqs. (2.11) and (2.13) in the Majorana 
neutrino case then follows from the relations —75^ = vl and —75^ = — vr. 

Consider now eq. (2.13). Using, as before, the block-matrix form for the neutrino 
propagator in the mixed coordinate-energy representation S(E;x,x f ) , we arrive at the 
equation 

—M* E + ia- V + V°- V-<x\ (S LL S LR \ = 3 _ , ft ()\ 

E-ia-V -V° -V-a -M J [s RL S RR ) l * X ' [o l] ' 

(A22) 

It differs from eq. (A2) by the presence of the potential-dependent term in the 12-entry 
of the first matrix on the left hand side. Since in the case of Majorana neutrinos M is in 
general complex symmetric, we replaced by M*. From (A22) we obtain a system of 
two coupled equations for Slr and Srr: 

-M*S LR (E; x, x') + (E + ia ■ V + V°(x) - V(x) • a)S RR (E; x, x') = , (A23) 
[{E - to- . V) - V°(x) - V(x) • a]S LR (E;x,x') - MS RR (E;x,x') = 5 3 (x - x') . (A24) 

Note that eq. (A24) coincides with (A4), whereas eq. (A23) differs from (A3) by an extra 
potential-dependent term in the coefficient of Srr(E;x,x.'). Next, we define, as before, 
J(E;x,x') = (M*)~ l Srr(E; x, x'). Eq. (A23) then gives 

S L r(E; x, x') = [E + ia ■ V + V°(x) - V(x) • tr] J(E; x, x') . (A25) 

For relativistic neutrinos (E+icr-'V)J(E; x, x') ~ 2EJ(E; x, x'), therefore, under condition 
(A12) one can neglect the term V°(x) — V(x) • a in (A25). 12 Eq. (A25) then reduces to 
eq. (A5). Since eqs. (A24) and (A4) coincide, we find that the neutrino propagator is still 
given by (A19), where F(E,x,x') satisfies eq. (A20) with the boundary condition (A10). 

Thus, for propagation of relativistic neutrinos in matter with potential satisfying 
|V^(x)| <C \E\ the propagator of Majorana neutrinos coincides with that of Dirac neu- 
trinos. 

12 Note that we cannot neglect the similar term in eq. (A24) because the coefficient of Slr(E;x,x') in 
this equation contains (E — icr ■ V) rather than (E + icr ■ V). 
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Appendix B: Matter-induced neutrino potentials 



We summarize here the expressions for the potentials of relativistic neutrinos caused by 
coherent forward scattering of neutrinos on background particles. For definiteness, we 
concentrate on the Dirac neutrino case; the potentials for Majorana neutrinos are the 
same, the only difference being that what we call antineutrinos in the Dirac case are just 
right-handed neutrino components in the Majorana case. 

Neutrino interact with matter through the charged current (CC) and neutral current 
(NC) interactions mediated by and Z° bosons, respectively. As we shall show, the 
effective Lagrangian of neutrino interaction with matter can be written as 

L mt = -v{ nt y^)P LV , (Bl) 

where the matrix of matter-induced neutrino potentials V 11 is diagonal in the flavour basis 13 
and is the sum of the CC and NC contributions: V 1 = Vq C + V^ c . Adding Ci nt to the free 
neutrino Lagrangian and making use of the standard Euler-Lagrange formalism to derive 
the neutrino equation of motion, one arrives at eq. (2.13) for the neutrino propagator in 
matter. 

We shall now concentrate on the potentials Vq C and V^ c . We will be assuming (ex- 
cept in eq. (B16) below) that the energies of neutrinos and particles of the medium are 
small compared to the VF-boson mass mw- I n an ordinary matter with no muons or 
tauons present, only electron neutrinos experience CC interactions, which are due to their 
scattering on the electrons of the medium. The effective Lagrangian of this interaction is 

£CC = -^f^)^ 1 " 75)e(x)][e(x) 7M (l - 75^0*0] , (B2) 

where Gf is the Fermi constant. We then employ the Fierz transformation to permute the 
neutrino field with the electron one and take the expectation value of the electron current 
over the state of the medium. This gives 

[CCCU = -VeJ» [(Ve)^ C (x)] P L V e , (B3) 

where [1, 16] 

(V e r cc (x) = V2G F (e(x)^(l - lb )e(x)) . (B4) 

Here (...) means the average over the state of the medium, and we have taken into account 
that for relativistic left-handed neutrinos (1 — 7s)^l ~ 2^l. By making use of the solutions 
of the Dirac equation for electrons, for the expectation values of the components of the 
electron current we find 

(e(x) 7 °e(x)) = N e (x) , (e(x)Ye(x)) = N e (x)vl(x) , (e(x) 7 °7 5 e(x)) = N e (x)(a e -v e ) , 

(e(z) 7 Se(x)> = JVe(a0[m e (ff7£ e ) + {[E e /(E e + m e )K(ovv e ))] , (B5) 

where N e (x) is the electron number density, v l e (x) is the ith component of the electron 
velocity, and a\ are the electron Pauli matrices (i = 1,2,3). Note that the expectation 



Except in media containing neutrino backgrounds, see below. 
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values of all the components of the axial- vector current vanish in a medium with unpolarized 
electrons. For such a medium from (B4) and (B5) we obtain 

(V e )° cc (x) = V2G F N e (x) , (V e y cc (x) = V2G F N e (x)vi(x) , (B6) 

The CC contribution to the expression V = V° — v„-V that enters into eqs. (2.27) and 
(A20) is then 

(V e ) CC = V2G F N e (l-v eC os6 eu ), (B7) 

where 9 eu is the angle between the momenta of the electron and the neutrino. For media 
with electrons at rest or non-relativistic electrons (v e <C 1) the spatial components of the 
CC potential can be neglected, and one obtains V C c - Vg c = y/2G F N e (x). This is the 
expression for the neutrino potential which is relevant e.g. for neutrino oscillations in the 
sun and inside the earth. It should be noted, however, that during supernova collapse or in 
rotating neutron stars bulk matter velocities may be substantial, leading to non-negligible 
net fluxes. In those cases the terms in the neutrino potentials that depend on the velocities 
of background particles should be retained. 

Consider now NC contributions to the matter-induced neutrino potentials. The effec- 
tive Lagrangian of the NC interaction of v a (a = e, /i, r) with a fermion / where f = e, p, n 
(or a background neutrino which may be abundant in supernovae or in the early universe) 
is 

£ NC = [vaW-fO- ~ 75K(x)] [ii f {x) lfM {T ZLf - 2Q f sin 2 9 w )^ f (x)] . (B8) 

Here Qf and T 3 £j are the electric charge of the fermion / and the third isospin projection 
of its left-handed component, respectively, and 9w is the Weinberg angle. Similarly to 
eq. (B3), upon averaging the variables of the fermion / over the state of the matter we find 

[£nc]v a ,f = -^alfj,[(V a )^ c (x)] f P L v a , (B9) 

where 

[(KOngO*)]/ = ^G F (T 3Lf - 2Q f sm 2 e w ) (^ / (x) 7 ^/(^)) • (BIO) 

It is important to note that the NC-induced potentials (BIO) do not depend on the neutrino 
flavour index a. That is, they are the same for all three active neutrino species (u e , Un 
and v T ) and vanish for sterile neutrinos. The equality [(V^jtc]/ = [(^m)nc]/ = [(^~)nc]/ 
actually holds only at tree level; at one-loop level tiny differences between these poten- 
tials arise, which are usually irrelevant. They can, however, play some role at extremely 
high densities, e.g. in supernovae. We will consider loop-induced NC contributions to the 
potentials below. 

The expectation value (V , /7^V'/) can be obtained from eq. (B5) by replacing the sub- 
script e by /. We will assume now that the particles / are unpolarized and have zero mean 
velocities or are non-relativistic and therefore one can keep only the time components of 
the NC-induced neutrino potentials. For the NC contributions of neutrino scattering on 



- 29 - 



the electrons, protons and neutrons of the matter we then find [17] 



[(V a hc(x)]e = V2G F N e (x)(-- + 2 S m 2 9 w ), 
[{V a )nc(x)] P = V2G F N p (x){^ - 2 sin 2 9 W 
l(V a ) NC (x)]n = V2G F (- (Bll) 

In an electrically neutral matter one has N e {x) = —N p (x), so that the electron and proton 
contributions cancel each other, and the only non-zero net effect is due to the neutrons. 
Combining the CC and NC contributions to the neutrino potential, in the flavour basis 
(z/ e , Vu, v T , v s ) where v s is a hypothetical sterile neutrino, we get for the matrix V = 
Vcc + ^NC 

V = V2G F dizg(N e -^, 0). (B12) 

For antineutrinos the right-hand side of this equality should be multiplied by —1. 

Since one can always add to the effective Hamiltonian in eq. (1.3) any matrix propor- 
tional to the unit matrix without affecting the oscillation probabilities, one can modify the 
matrix V in (B12) according to V ->■ V + G F (N n /y/2)-l. This yields 

V = V2G F diag(N e , 0, 0, (B13) 

For neutrino propagation in normal media, this form is the most often used one. It is 
especially convenient when oscillations between only active neutrino species are considered, 
since in this case the matrix V in (B13) has only one non-zero element V e (x) = \^2G F N e (x). 
Thus, in this case only the CC contribution to the neutrino potential affects neutrino 
oscillations. Note a useful relation \^2G F N e ~ 7.63 x W~ 1A pY e eV, where p is the matter 
density in g/cm 3 and Y e is the electron fraction (number of electrons per baryon) in matter. 

Loop corrections to the matter-induced neutrino potentials were calculated in [18]. 
They differ for neutrinos of different flavour due to the differences of the masses of the 
corresponding charged leptons. The most important difference is the one between the 
potentials of v T and u^, since this difference vanishes at tree level (see (B12) or (B13)). For 
a neutral unpolarized medium it is 



3 



m 2 w ,„ 



AK, = Vr-Vfi~ ±^G F m 2 T (N p + N n ) In -f - (N p + -N n ) . (B14) 



ml v F 3 



Here and below the upper sign always refers to neutrinos and the lower one to antineutrinos. 
Note that AV Tfl is very small, AV T u/V e ~ 5 x 10 -5 . However, it may play some role at 
very high densities, in particular, for supernova neutrinos [19]. One-loop contributions to 
AVt^ in the neutrino backgrounds were calculated in [20]. 

The above formulas for the matter-induced neutrino potentials apply to the ordinary 
unpolarized matter at zero temperature and with no antiparticles. We will now relax these 
constraints. 
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Neutrino potentials in hot and dense matter and in neutrino backgrounds 

This case is relevant for the early universe and supernova physics. It was studied in 
refs. [17, 21], the results of which we summarize here. In an electrically neutral unpolarized 
medium consisting in general of electrons, muons, r-leptons, protons, neutrons and their 
antiparticles with zero mean velocities, the potential of electron (anti)neutrinos is 

r *2E 
V e = ±V2G F (N e -N s )--(N n -N n )T — ((E e (l+v 2 e /3))N e + (E iS (l+vl/3))N s ) . (B15) 

Here E is the energy of the neutrino, E e and Eg are those of the electrons and positrons 
of the medium, v e and Vg are the electron and positron velocities, and iVj stands for the 
number density of the antiparticles of /. All the averages are now taken over the proper 
thermal distributions of the background particles. Note that the first and second terms in 
the square brackets in (B24) are the generalizations of the CC and NC contributions to 
V e discussed above to the case when antiparticles are present in matter. The NC-induced 
term comes only from the neutrino scattering on neutrons, since the NC contributions of 
all charged particles cancel in an electrically neutral medium. The third term in (B24) is 
due to CC and is rather special. It comes from the second-order term in the expansion of 
the propagator in powers of l/mfy. Due to an extra power of in the denominator, 
it is negligibly small in an ordinary matter. However, it does not vanish in the limit 
iV e = iVg and so becomes important in a medium with equal (or almost equal) abundances 
of particles and antiparticles, when the contributions of the first two terms are negligible. 
In addition, this term has the same (negative) sign for electron neutrinos and antineutrinos. 

The last property, as well as the fact that the third term in (B15) is non-zero for 
iV e = Ng, can be understood as follows. The contribution of the W-boson exchange to 
v e e scattering amplitude is proportional to g 2 j\m^ v — (q — p) 2 ] « (5 2 /m^ / )(l — 2q-p/m 2 v ), 
where q and p are 4-momenta of the neutrino and of a background electron, and g is the 
SU(2)l gauge coupling constant. For v e scattering on positrons one has to flip the overall 
sign of this expression and to replace p — > —p, so that the corresponding contribution to 

V e is proportional to — g 2 /[ m w/ ~~ (.1 + P) 2 \ ~ ~(d 2 / m w)(.^ + ^Q'P/ m w)- Obviously, the 
terms ~ l/rriy^ enter with the same (negative) sign. The situation is similar if one goes 
from neutrinos to antineutrinos, in which case one has to replace g 2 — > —g 2 , q — > —q. The 
factor (E e (l + v 2 /3)) comes from the averaging of E e (l — v e cos#q P ) 2 over the angle # qp 
between the momenta of the neutrino and the background electron. 14 

Another interesting propagator effect takes place at extremely high neutrino and/or 
electron energies. In a CP-symmetric matter with equal electron and positron abundances 
the CC contribution to the matter-induced self-energy of v e is proportional to 



9 — 77 T?, (B16) 



14 One power of (1 — v e cos# qp ) is due to the fact that for relativistic neutrinos V = V — V 3 ~ V°(l 
v e cos#q P ), while the other power and the factor E e come from p-q ~ EE e (l — v e cos^qp). 
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where no expansion in powers of l/m^y has been done. In the limit (p-q) 2 <C the 
previous results are recovered, whereas we see that for 4(p-g) 2 > the potential changes 
its sign. 

The potentials of v^ipu) and u T (p T ) in matter are given by expressions similar to (B15), 
with the index e replaced by [i or r, respectively. If the medium contains no fft 1 and r ± , 
only the neutron contributions to and V T (which coincide with the second term in (B15)) 
survive. 

In a number of applications (e.g. for supernova neutrinos) it is necessary to consider 
neutrino potentials in neutrino backgrounds. Those are due to the NC interactions, and 
they depend on whether the background of the same flavour or different flavour neutrinos 
is considered. In the case of the same flavour neutrino background the corresponding 
contribution to the potential of the test neutrino of momentum q is 

AV a = V2G F J4^{ ± 2«(p) - <(p))(l - cos qp ) 

- ^# k(p)<(p) + E Pa (p)4 a (p)] (1 - cose qp ) 2 }. (B17) 

Here n^ a (p) and rip (p) are the occupation numbers of the left-handed background neutri- 
nos of flavour a and of their antiparticles. The quantity (p) is related to the neutrino 
number density through 



(ftp L 

and similarly for antineutrinos. The origin of the last term in the curly brackets is similar to 
that of the last term in (B15), except that it comes from the expansion of the Z boson rather 
than W boson propagator. Note that the neutrino potential due to the coherent forward 
scattering on a background neutrino vanishes when the velocities of the test and background 
neutrinos are parallel to each other, i.e. when cos# qp = 1. This happens because there is 
no forward neutrino-neutrino scattering for completely relativistic neutrinos moving in the 
same direction. If the momentum distribution of the background neutrinos is isotropic, 
then (cos#q P ) = and (B17) reduces to 



AV a = ±2V2G F (Ni - At) - 8E »°Wf G r 

z 



(E Ua )Nt + (E Pa )Nl(p) . (B19) 



For a test neutrino in a neutrino background of different flavour one has 
AV a = ±V2G F J ^3 (n L V!i (p) - n L Ufj (p)) (1 - cos qp ) (/?/«). (B20) 

The extra factor of 2 in front of the first term in (B17) in comparison with (B20) is 
due to the exchange effects in the case of same-flavour neutrino background. If the 
momentum distribution of the background neutrinos is isotropic, eq. (B20) reduces to 
AV a = ±V2G F (N^-N^). 
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Unlike in ordinary matter, neutrino potentials in neutrino backgrounds are not in 
general diagonal in the flavour basis. While the diagonal terms (B17) and (B20) arise 
from the coherent forward scattering processes f a (k) + z^(p) — > f a (k) + vp{p) (where the 
neutrino momenta are shown in the parentheses), the NC-induced momentum-exchange 
processes v a (k) + vp{p) — > f«(p) + ^(k) with a ^ (3 are also coherent and lead to flavour- 
off-diagonal potentials V a p [21-24]. The potential V a p due to the scattering of a test 
neutrino of momentum q on background neutrinos and antineutrinos is 

V a p = V2G F J-0-{(p^(p) - p^»)(l - cos V) 

(l-cos0 qp ) 2 }- (B21) 

Here p„ al//3 {p) and Pp a p p {p) are the off-diagonal elements of the density matrices of left- 
handed neutrinos and their antiparticles in the flavour space: 

p^(p) = <4l(pK^(p)>> ( B22 ) 

where a^iip) and a a L(p) are the production and annihilation operators of iy a i(p), and sim- 
ilarly for antineutrinos. Note that the neutrino occupation numbers that enter in eqs. (B17), 
(B18) and (B20) are the diagonal elements of these density matrices: n^ a (p) = p^ aPa {p), 
n P a (p) = Pp a D a (p)- The off-diagonal potentials V a p are in general complex, with Vp a = V*a. 
Eq. (B21) is valid for test neutrinos; for antineutrinos one has to replace 

Pu a u,iP) ^ pLp,(p) , E Va {p) -+ E Pa (p) , ^(p)«^(p). (B23) 

When considering neutrino flavour evolution in matter, one usually assumes that there 
is no back reaction of this evolution on the properties of the medium, and therefore matter- 
induced neutrino potentials are fixed external quantities. This is in general not true for 
neutrino oscillations in neutrino backgrounds, as the oscillations affect the state of the 
background. Therefore describing neutrino oscillations in media containing significant 
abundances of background neutrinos represents a complex non-linear problem. The el- 
ements of the neutrino and antineutrino density matrices in the flavour space that enter 
into eqs. (B17), (B20) and (B21) must then be found self-consistently as solutions of the 
same flavour evolution problem. 

Magnetized matter 

In a medium with a magnetic field the particles of matter have in general non-zero average 
spin. In this case one can no longer neglect the axial- vector contributions to the neutrino 
potentials (see eq. (B5)). Under realistic conditions the average spin of the particles is 
relatively small, so that their polarizations are linear in the magnetic field strength. In this 
case in a matter consisting of electrons, protons and neutrons the neutrino potentials V® 
get the extra contributions 

AV£=±(c<i v + c% + tf> + c n )B ]] , 

AVl tl/T = ±(c* z + c? + c n )B ll , (B24) 



EvAp)pI 



i I'd 



{p) + E Vfj {p)p^ aUfj {p) 



- 33 - 



where Bu is the component of the magnetic field along the neutrino velocity. The coef- 
ficients Cyy and c z describe the contributions to the neutrino potentials coming from the 
polarization of the background electrons and caused by the CC and NC interactions re- 
spectively. The coefficients c p and c n are due to the polarization of the background protons 
and neutrons. For a relativistic gas of degenerate electrons (i.e., for Ep 3> T where Ep is 
the electron Fermi energy and T is the temperature), such as e.g. in or near the supernova 
core, one has [25-27] 

4-?%(^y /3 , ^ = -2c|. (B25) 



2^/2 

For the contributions of the polarization of non-relativistic protons and neutrons with 
Boltzmann distributions functions one finds [28, 29] 

cP ^ ^&M!L Np , e . ^ 9 l^N n . (B26) 

Here /xjv = e/(2m p ) — 3.152 x 10~ 18 MeV/G is the nuclear Bohr magneton, fj, p and n n 
are the proton and nucleon magnetic moments in units of the nuclear Bohr magneton 
(fi p = 2.793, /i n = —1.913), and g v A and g A are the NC axial-vector coupling constants 
of proton and neutron. For free nucleons, one has g A ~ 1.36 and g A ~ —1.18 [31, 32]. 
In applications for neutron stars, the values of g^ A and g A in nuclear matter are more 
relevant; they can be estimated as free-space values divided by 1.27 [31], i.e. g p A « 1.07, 
g A ~ —0.93. Note that c p , c n and c z are all of the same sign. For non-degenerate particles, 
thermal fluctuations tend to destroy the polarization, and therefore cP and c n decrease with 
increasing temperature T. 



Appendix C: Proof of the equalities p* = p^- and p'^ = p/</ 

We shall prove here that for macroscopic distances |x' — x| the momentum integrals in the 
expression 

V(£;x',x) = J ^ Sp a (E; p', p) j^-ip* , (C1) 

receive their main contributions from small regions around the points p = p* and p' = p*, 
which are defined as follows. For a given E, the value of p* is obtained from the dispersion 
relation that stems from the neutrino evolution equation in matter of constant density 
equal to the density at the initial point of neutrino evolution x. Likewise, p'^ is found from 
the neutrino dispersion relation in matter of constant density corresponding to the final 
point of neutrino evolution x'. 

Let us first consider Sp a {E\ p', p), which is a Fourier transform of Sp a (E; x', x) (see 
eq. (2.24)): 

Sp a (E;p',p) = J d 3 xd 3 x'Sp a (E;x',x)e- ip ' x ' +ipx 

= -^f d 3 xd 3 x / ^^F^( J E;x / ,x)e- ip ' x ' +ipx+i l i? ll x - x 'l . (G2) 
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Here in the second line we used eq. (2.26). For typical values of the energy E and momenta 
p and p' of interest to us, the integrand in (C2) contains a fast oscillating phase factor, 
and therefore the integral can be calculated in the stationary phase approximation (see, 
e.g., [30]). Defining 

G(x',x) = -p'x' + px+ |£||x-x'| - ilnFp a (E;x',x) , (C3) 

we can rewrite eq. (C2) as 

%*(£;p',p) = d 3 xd 3 x'-^-e iG ^'^. (G4) 

The main contributions to the integrals over the coordinates come from small neighbour- 
hoods of the points where the phase G(x',x) is stationary. These points are found from 
the conditions 

VG(x,x) = -p +i\E\r-i , =0, (C5) 

Fp a (E;x',x) 

VG(x,x) = p-\E\r-i ^ = ' ( C6 ) 

Fp a (E;x,x r ) 

where V' is the gradient with respect to the coordinate x' and f = (x' — x)/|x' — x|. 
Eqs. (C5) and (C6) can also be rewritten as 

iV'Fp a (E;x',x) = -(p'- \E\v)F pa {E-x',x) , (C7) 
iVFp a (E;x',x) = (p-\E\r)Fp a (E;x',x). (C8) 

Eqs. (C7), (C8) (or (C5), (C6)) should be solved with respect to the coordinates x and 
x' for fixed values of p and p'. We denote the corresponding solutions x* and x^. Notice 
that x* and x' # , are functions of p and p'; we will not indicate this dependence explicitly 
in most of the following formulas in order not to overload the notation. 
Applying the stationary phase approximation to eq. (C4) yields 



Sp a (E; P >,p) « & E ,, J , n ff M e^<^ , (C9) 

2tt|x* - x;| y |L)(x;,x*)| 

where rj is a constant phase which is of no relevance for us, and 

Z)(x^,x*) = det 



r/9 2 G(x',x)\ 




A dx^dxj J 





(CIO) 

Next, we substitute (C9) into (CI). Since for macroscopically separated x and x' the 
integrand of (CI) contains a fast oscillating phase factor (see section 2.2), we can calculate 
the integrals over p and p' by once again making use of the stationary phase approximation. 
In doing so, we will need to find the stationary points of the expression 

G(p / ,p)^ P V-px + G(x;,x,). (Cll) 
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Here we have taken into account that Z)(x^,x*) is not a fast oscillating function and 
therefore, in keeping with the stationary phase approximation, it need not be included in 
the phase factor G(p, p') but can instead be left as a pre-exponential factor. Substituting 
(C3) into (Cll) yields 



G(p',p) = p'(x' - x'J - p(x - x*) + |£||x* - x'J -iln#fl a (E;xl,x*) 



(C12) 



Let us now find stationary points of G(p',p), which will give us the momenta that 
yield dominant contributions to the integrals over p and p' in (CI). Requiring that the 
derivatives of G(p', p) with respect to the components of p vanish, we find 



dG(p',p) 

dpi 



[OC X* 



dx 



dpi 



+ 



Fp a (E;x'^x 
d 







dx 



*J 



' dx 



— i^Q,(E;xi,x*) + (p'j - \E\rj)Fp a (E;x.'*,x.*) 



*.i 



dx' 



dpi 



(C13) 



From eqs. (C7) and (C8) it follows that the expressions in square brackets in (C13) vanish, 
so that (C13) simply reduces to x = x*(p, p'). Quite analogously, by requiring that the 
derivatives of G(p', p) with respect to the components of p' vanish, one finds x' = x^(p, p'). 
Thus, the momenta at which the phase G(p', p) is stationary are obtained as the solutions 
of the system of equations 



x' = x*(P,p'), 



x = x*(p,p 



(C14) 



We will call the corresponding solutions p* and p^. Recall now that x* and x^ are the 
solutions of the system of equations (C7) and (C8) for fixed values of p and p'. From 
eq. (C14) it follows that p* and p^ are the solutions of the same system (C7), (C8) which 
should now be considered as equations for the momenta at fixed values of the coordinates 
x and x'. Note that when considered as equations for the momenta, eqs. (C7) and (C8) are 
actually much simpler than when considered as equations for the coordinates; for known 
F/3 a (E; x', x) one finds the solutions for the momenta p and p' immediately - they are 
simply given by (C5) and (C6). 

Next, we recall that the components of p* and p'* that are orthogonal to the vector 
x' — x are negligibly small in all situations of practical interest (see discussion in section 2.2); 
therefore we are only interested in the longitudinal components of these momenta, which 
we denote p* and p'*. Multiplying eqs. (C5) and (C6) by rFp a (E;x,x') yields 

i-^F 0a {E;x' ,x) = -( P ; - \E\)Fp a {E;x',x) , 



i—F 0a {E;x',x) = (p* - \E\)Fp a {E;x',x), 



(C15) 
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where d/dx' = f -V' and d/dx = f • V. On the other hand, we have 

i^jHE; x', x) = ff(x') F(E; x', x) , (C16) 

i-^F(E; x', x) = —F(E; x', x) H(x) . (C17) 

ax 

where -ff(x) = MM'l'/2|£'|+y(x). The first of these equations is just eq. (2.27), whereas the 
second one, which involves the differentiation of F(E;x' ,x) with respect to the coordinate 
of the initial rather than final point of neutrino propagation, can be derived from the first 
one. 15 Using eqs. (C16) and (C17) in (C15), we find 

(pi - \E\) F(E; x', x) = -H(x>)F(E; x', x) (C18) 
(p* - \E\) F(E; x', x) = -F(E; x', x)iT(x) . (C19) 

Since the effective Hamiltonian H (x) is non-diagonal in the flavour-eigenstate basis, eqs. (C18) 
and (C19) are matrix equations for p* and p'. They are simplified in the local matter eigen- 
state bases defined in eqs. (2.29) and (2.30). In these bases the effective Hamiltonians H. at 
the initial and final points of neutrino propagation are diagonal: H(x.)km = H{x)k8km, 
H(x.')k'M' = 7-L{~x!)k>Sk'M>- Here Hk^) is the Kth. local eigenvalue of H at the point 
with the coordinate z. Thus, we finally obtain from (C18) and (C19) 

p'*=v' K , = \E\-n K > (*■'), (C20) 
p* = VK = \E\ - H K (x) ■ (C21) 

Eqs. (C20) and (C21) give the longitudinal (with respect to x' — x) components of the 
vectors p* and p^; as discussed above, their transverse components nearly vanish: 

P*±-0, p'^O. (C22) 

Note that eqs. (C20)-(C22) yield the correct neutrino dispersion relations in the limits of 
vanishing vacuum mixing or vanishing matter density. 

Thus, we have proved that main contributions to the momentum integrals in (CI) 
come from small regions around the of momenta p* and p^, which satisfy the dispersion 
relations in matter at the initial and final points of neutrino propagation, respectively. 

Appendix D: Evolution equation in the adiabatic regime 

We shall prove here that in the adiabatic regime, when matter density varies sufficiently 
slowly along the neutrino path, the amplitude of the overall neutrino production-propagation- 
detection process (4.6) satisfies the standard evolution equation (1.3). 

15 Indeed, F(E;x' ,x) can be written as F(_E;x',x) = F(E; x', xi)F(E; xi, x) with arbitrary xi. This 
relation can be easily verified by substituting it into eq. (C16). Then, from (d/dxi)F(E; x', x) = we have 
[(d/da:i)F(£;x',xi)]F(£;xi,x) = -F(E; x', xi)(d/da;i)F(£; xi, x). Substituting here (d/dxi)F(E; xi, x) 
from (C16) and multiplying the result by [F(E; xi, x)] _1 on the right, one arrives at (C17). 
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In the adiabatic regime the transitions between different matter eigenstates are strongly 
suppressed, i.e. all matter eigenstates evolve independently. This means that the quan- 
tity T that characterizes neutrino propagation in the matter eigenstate basis is diagonal: 
Fk'k{E; x, Xo) = Tk(E; x, Xo)Sk'k- The amplitude (4.6) can then be written as 

^(£,x,x ) = {C/(x)[^;x,x )<I>P^]C/ t (xo)} / 3 Q , (Dl) 
where all the factors in the square brackets are diagonal. From eq. (2.32) we have 

f(E; x, x ) = tf(x)F(E; x, x )L>(x ) , (D2) 
so that (D2) can be rewritten as 

A%(E,x,xo) = {F(i?;x,xo)L>(xo)«I>P^C/ t (xo)} /3Q , (D3) 
Substituting this into eq. (Dl) and differentiating, we obtain 

i^ tot (£, x, x ) =i-^{F(E; x, xo)L>(x )$P^C/(xo) t } 

=H(x)F(E; x, xo)f/(xo)«»p$ D [/(xo) t = H(x)A tot (E, x, x ) , (D4) 

where we used eq. (2.27). 

Thus, in the adiabatic regime the amplitude of the overall process satisfies the standard 
evolution equation (1.3), irrespectively of whether or not the conditions of coherent neutrino 
production and detection are satisfied. However, the boundary condition for this amplitude 
differs from the standard one. Instead, from eqs. (D3) and (A10) we find 

^£(£,x,x )U X0 = Y J U(^T aK Up K {^p(E, VK )^ D (E, V ' K ) = {U^pQdU^xo)} 

K 

(D5) 

If the coherence conditions for neutrino production and detection (3.12), (3.13) are satisfied, 
one can replace the momenta p^ and p' K in the arguments of the amplitudes <3>p and 
3>£) in eq. (D5) by the corresponding average values and pull these amplitudes from the 
sum. Eq. (D5) then reduces, up to a constant factor, to the standard boundary condition: 
x.xoJIx-mco = 8p a Qp(E,v)* D (E,tf). 
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